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Abstract 

Using perturbative methods we derive new results for the spectrum and correlation func- 
tions of the general 23-chiral Potts quantum chain in the massive low-temperature phase. 
Explicit calculations of the ground state energy and the first excitations in the zero mo- 
mentum sector give excellent approximations and confirm the general statement that the 
spectrum in the low-temperature phase of general 2 n -spin quantum chains is identical to 
one in the high-temperature phase where the role of charge and boundary conditions are 
interchanged. Using a perturbative expansion of the ground state for the Z3 model we are 
able to gain some insight in correlation functions. We argue that they might be oscillating 
and give estimates for the oscillation length as well as the correlation length. 
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1. Introduction 



The self-dual Z3-chiral Potts model was introduced by Howes, Kadanoff and den Nijs [1] 
and studied using e.g. fermionization and approximative methods, in particular pertur- 
bation expansions. One remarkable result of the perturbation expansions was that for 
special values of the parameters the first translationally invariant excitation is linear in 
the inverse temperature A for special values of the parameters, v. Gehlen and Rittenberg 
then realized [2] that this model is integrable for these special values of the parameters 
because it satisfies the Dolan-Grady integrability condition [3] which is equivalent [4] [5] to 
Onsager's algebra [6]. They also generalized this 'superintegrable' chiral Potts model to 
general 2 n -spin n [2] . Afterwards, it attracted much attention because it can be related to 
a classical model that satisfies a generalized Yang-Baxter equation with Boltzmann weights 
defined on higher genus Riemannian surfaces [7 — 11]. However, even then perturbative 
methods lead to important new results [12] [13]. One example is a conjecture for the exact 
form of the order parameters in general superintegrable Z n -chiral Potts chains [13]. 

Recently, a particle interpretation of the momentum zero sectors in the high-temperature 
phase of all 2 n -chiral Potts models at general values of the parameters has been proposed 
[14] and a quasi-particle spectrum has been derived for the superintegrable 23-chiral Potts 
model [15]. Furthermore, a scaling exponent for the wave vector in the low-temperature 
phase of the 2 3 -chiral Potts model has been calculated in [16] [17] from level crossings in 
the ground state. This motivated us to perform the perturbative calculations reported 
in this paper. On the one hand the excitation spectrum in the low-temperature phase 
of general 2 n -chiral Potts quantum chains is not completely understood. On the other 
hand very little is known about correlation functions. Perturbation expansions enable us 
to gain more insight in the spectrum of the low-temperature phase and shed some light on 
correlation functions. 

A general 2 n -spin quantum chain with N sites is defined by the Hamiltonian: 

N n-1 

h n ] = - E E ~ x ^ a J + ^ r J r l+i- (i-i) 

j = l k=l 

For n = 2 (1.1) is just the well known Ising model. In this case the operators <jj and Tj 
are the Pauli spin matrices a x and a z acting in a vector space (T. 2 located at site j. For 
general n one may think of the operators Oj and Tj as generalizations of the Pauli spin 
matrices - see (2.3) below. In this paper we identify the iV + 1st site with the 1st site, i.e. 
we use toroidal boundary conditions. 

The Hamiltonians (1.1) contain 2n — 1 parameters. The temperature like parameter 
A will be chosen real while the coupling constants a k and a k will be generally complex, 
is hermitian iff a k = ot* n _ k and a k = OL^_ k . 

In this paper we will parametrize the constants a k , a k in (1.1) by two angles 0, cp, fixing 
their dependence on k: 

e ^-i) e ^(f-i) 

«fc = > a k = . „ k ■ (1-2) 

sin — sin — 

n n 
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(1.2) is called the general 'chiral Potts model'. This parametrization is convenient because 
it can easily be specialized to various models. Setting = <p = yields models with a 
second order phase transition at A = 1 that can be described by a parafermionic conformal 
field theory in the limit N — > oo [18]. These so-called Fateev-Zamolodchikov-models lead 
to extended conformal algebras WAn-i where the simple fields have conformal dimension 
2, . . . ,n [19] [20]. The spectrum of the Hamiltonian (1.1) with = (p = can be described 
by the first unitary minimal model of the algebra WA n -i- For n = 3 it coincides with 
the three-states Potts model and the symmetry algebra is Zamolodchikov's well known 
spin-three extended conformal algebra [21] at c = |. 

For = (p = ^ (1.2) specializes to the 'superintegrable' 2 n -chiral Potts model which 
exhibits remarkable integrability properties [2]: At = V? = § the Hamiltonians (1.1) 
satisfy the Dolan-Grady integrability condition [3]. 

Albertini et al. [7 — 9] have shown that the Hamiltonian (1.1) can be obtained for more gen- 
eral values of the angles 0, <p as the r-continuum limit of an integrable classical statistical 
model if one imposes the constraint 



H^' with the choices (1.2), (1.3) is in general no more self-dual. However, if we choose 

= (p in (1.2) is self-dual. Sometimes (1.3) is implied when referring to the chiral 

Potts model but we prefer to call (1.1) with (1.2) the general chiral Potts model. 

The Z3 version of (1.1) is known to have four phases [8] [22] [17]: Two massive and 
two massless phases. One of the massive phases is ordered and the other massive phase is 
disordered. The low-temperature phase (small A) that we are going to study in this paper 
is the ordered massive phase. At = ip = ^ it appears in the range < A < 0.901292 . . . 



In the next section we will review some well known facts about the Hamiltonian 
(1.1). Then, in section 3 we will evaluate the ground state energy and the first excitations 
perturbatively for n = 3 and zero momentum, but arbitrary 0, (p. The results remind us 
of a general duality statement that we discuss in section 4. Finally, section 5 is devoted to 
a study of correlation functions of the Z 3 -chain. 

2. Preliminaries 



In this section we summarize well known facts about Z n -quantum spin chains and introduce 
notations that will be useful later on. For a more detailed, recent review see e.g. [23] . 

First, we give a precise definition of the operators Tj and aj. <jj and Tj freely generate a 
finite dimensional associative algebra by the following relations (1 < j,l < N): 



A cos ip = cos 0. 



(1.3) 



[8] [22]. 



TjTi = TiTj , 



n pn n 



(2.1) 
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where u is the nth root of unity uj = . In this paper we will only consider boundary 
conditions of type r_/v+i = u~ R Ti, R e Z n for the Hamiltonian (1.1). We will mainly 
focus on periodic boundary conditions T^ + i = Ti. 

The algebra (2.1) can be conveniently represented in ® N (L n . In this space we can 
choose the following basis if we label the standard basis of (L n by {eo, . . . , e n _i}: 

\ix . . .i N ) := e ix ® . . . <g> e iN 0<ij<n — l. (2.2) 

Usually one considers a special representation r of the algebra (2.1) - a definition can 
be found in appendix A. However, for low-temperature expansions of (1.1) it is more 
convenient to consider a different representation f: 

r(Tj) | ii . . . ij . . . i N ) = u l > \ h . . . ij . . . i N ) ^ 
f(o"j) | i\ . . . ij . . . in) =| ii ■ ■ ■ (ij — 1 mod n) . . .in) ■ 



commutes with the Z n charge operator Q = Ylf=i a h ^ nus nas n cnar g e sectors. 
The eigenvalues of Q have the form u Q with Q e {0, . . . , n — 1}. We will refer to the 
number Q as the 'charge'. 

also commutes with the translation operator T^. The eigenvalues of T/v are iVth 
roots of unity. We label them by e tP and call P the 'momentum'. For a chain of length 
iV one has P G {0, . . . , 2?r ^ 1 ' > }. The eigenstates ||zi . . . in)) p with momentum P can 
be obtained from | i\ . . .in) by finite Fourier transformation. Tn acts for R = on the 
states (2.2) as 

f (T/v) | iii 2 ...in) =| ii ■ ■ ■ iNh) ■ (2.4) 
In the case R = the eigenstates with momentum P are given by 

N-l 

— i N )) P :=^E e iPx f(T N )- x \h...i N ). (2.5) 

Af is a suitable normalization constant. If the state \ h ■ ■ - In) has no symmetry, one has 
Af = N. This will apply to all cases below where we need (2.5). For a definition of the 
momentum eigenstates (2.5) in case of boundary conditions R ^ see e.g. [24] [23] 

3. Spectrum in the low-temperature phase 

In this section we calculate the ground state energy and the lowest excitations in the low- 
temperature phase using perturbative expansions according to [25] around A = of the 
Hamiltonian (1.1). We restrict ourselves to periodic boundary conditions R = 0. 

In [1] perturbation series were presented for the disorder operator (or magnetization m) and 
the first energy gap in the momentum zero sector of the superintegrable Z3-chiral Potts 
model leading to exact conjectures for both of them. After the superintegrable chiral 
Potts model had been generalized to general n [2] perturbation series for the ground state 
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energy, energy gap in the momentum zero sector, magnetization and susceptibility of this 
superintegrable 2 n -chiral Potts model were presented in [12]. At the same time elaborate 
expansions of the ground state energy and some excitations of the superintegrable chiral 
Potts model for n G {3, 4, 5} and in particular perturbation series for the order parameters 
with general n were calculated in [13]. First perturbative results for the energy gaps at 
more general values of the angles 0, <p were obtained in [14] where second order high- 
temperature expansions for the translationally invariant energy gaps in each charge sector 
of the general self-dual Z3- and Z4-chiral Potts models as well as a first order expansion 
for the dispersion relations for general n was presented. 

In this section we restrict once again to the ^-version of the chiral Potts model (1.1) but 
impose no restrictions on the angles 0, ip. We present low-temperature expansions for the 
ground state energy and in particular the first translationally invariant energy gaps that 
up to now have not been treated by perturbative methods because of high degeneracies. 

The normalization of the Hamiltonian (1.1) is chosen such that expansions around 
zero temperature A = are possible. If one wants to calculate expansions around infinite 
temperature one usually normalizes the Hamiltonian = , sets A = A -1 and 

performs expansions around A = 0. 

In each charge sector Q of the low-temperature phase there is one unique ground state. 
For arbitrary n it is given by: 

n— 1 

|GS;Q> :=-=Y. J ' Q \ 1 '-- 1 ) 

v 71 i=o 

provided that — \ < < § • For n = 3 (3.1) is the ground state if — tt < < ix and for 
n = 4 (3.1) is the ground state for — ^ < < The excited states are more complicated 
and highly degenerate. The space of the first excitation is spanned by those states which 
have precisely two blocks of different spins. Furthermore, the values of the spins in these 
two blocks must have difference one. For fixed P, Q and — f < < f we can choose the 
following basis for the space of the first excitation: 



\ a k 



n — l 

Q ) := -= Vc/ Q || (Z + l mod n)...(l + l mod n)l.. .1)) P . (3.2) 



In order to perform explicit calculations we now specialize to n = 3 with P = 0. The 
constant contribution in A to the ground state energy E^q.q^ and the first gap AEq^i can 
be calculated easily: 

4%Q> = -»js » s (!) ■ A£ S = *^<™ (!) • (3-3) 

As far as the ground state is concerned, we notice from the explicit form of the potential 
in (1.1) that the first order in the perturbation expansion vanishes. The next order that 
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has to be calculated is the second order E^ s _^ . One applies the potential V = J2j k a k°~j 
to the ground state: 

f(V) \GS;Q) = ( e * | a?) + e-*J> . (3.4) 

Using AEq \ = E^q^ — Epg.g^ one obtains the result 



p (2) v IH WI GS ;Q) I 2 

^PS;Q> p (0) _ p (0) 

|a)/|GS;Q) ^PS;^) ^|a) 

| (a? | f (V) | GS; Q) | 2 + | (a%_ l \ r(V) | GS; Q) | 2 



(3.5) 



2iV 



Sv^cos (|) 

Set C := cos (3)5 ^ := 1 — 4C 2 . Then, calculating higher orders in the same manner, one 



arrives at: 



v^os^f 3 + 4 l X 5| +0(A6) _ 



81C 2 4C 2 ?e 



(3.6) 



for sufficiently large N. For </? = = ^ (3.6) reproduces the result of [12]. The orders eg 
of the free energy per site eo are independent of iV if N > k. This is a general feature of 
the ground state energy for spin quantum chains with nearest neighbour interaction. 

Note that the expansion in powers of A (3.6) of eo does not depend on the charge sector 
for large N. From the explicit calculations we see that this is a general result: The order 
X k e^ of the free energy does not depend on the charge if iV > k. However, ground 
state level crossings have been observed in [16] [17]. Indeed, for short chain length iV high 
orders X k (k > N) do depend on the charge Q. Thus, at fixed A in the massive phase 
level crossings in the ground state do occur although the Q-clependent term of the ground 
state energy e decreases fast in magnitude with increasing N. The presence of such level 
crosssings is a hint for oscillating correlation functions and, in fact, the critical exponent 
of the wave vector can be calculated from them [26] [16] [17]. For the superintegrable case 
<p = tp = 77 we explicitly determined level crossings in the ground state using 9th resp. 
10th order expansions for 3 < iV < 7 sites. Our results for the temperatures A where the 
gap vanishes are in good agreement with the values presented in [17] (Table 3) and show in 
particular that there are no crossings besides the one presented in [17]. However, we have 
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argued that for larger iV we would need even higher orders for the study of level crossings 
which goes beyond current computer power. Thus, numerical methods provide a much 
better tool for the study of level crossings [16] [17]. 

Finally, we should mention that the approximation (3.6) is excellent for the whole massive 
low-temperature phase up to its boundary near A = 1 already for moderately long chains. 
For example at the phase transition A = 1 of the parity conserving Potts case = ip = 
one observes only a deviation of 0.8% between (3.6) and the exact result for N = 12 sites. 
In the superintegrable case (p = p= ^a,t\ = 0.9 (close to the boundary of the phase) 
the deviation of (3.6) at N = 12 from the exact result is also as small as 1.3%. Thus, 
the perturbation expansion (3.6) yields a surprisingly good approximation in the complete 
massive low-temperature phase. Even at the boundary of the phase the deviation is smaller 
than 2% for N > 10. 

,2 

In principle, one can derive a critical exponent a for the specific heat from a per- 
turbation expansion of the ground state energy cq. However, (3.6) leads only to a third 
order expansion for the specific heat and on the superintegrable line even the third order 
vanishes. Although the approximation of (3.6) to the ground state energy eo itself is so 
good, one certainly needs higher orders for accurate estimates of the critical exponent of 
the specific heat. Consequently, a has been estimated using a 13th order expansion of eo 
in [1] for the self-dual case <p = <p = and for the superintegrable case = tp = ^ the 
Ising-like form of the eigenvalues has been exploited to calculate even higher orders of e 
in [13]. The results in [1][13] indicate a = | independent of the angles 0, (p. 

The calculation of the smallest gap AEq^i is more difficult. Let q be the projector 
onto the space spanned by the states | a®) (3.2). In this space, the potential acts as follows: 



qr{V) \ a?) 


2 iV 

-7i (2eT| 


a?) +e*a;«|ag_ 1 >) 




qr{V) |a?> 


2 i<P . 

= "7! (2eT 1 


a? +1 > +2e"^ lag.!)) K k < N - 1 


(3.7) 


qr(V) |a^_!> 


= -7I (2e ~" 


l4_ 2 >+e-^|a?>). 





In the limit iV — > oo the lowest perturbative eigenvector converges to 

1 N-l 

Using (3.7) and (3.8) we can calculate that for iV — > oo and < f 

^lim AE Q , 1= 4V3C-A-|cos(|) 

p [ 8C(l + cos(f)) 4 (2 - cos (f)) | ~ 3 
\ 3V3(TZ + 2) 3V3C J 
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(3.8) 



(3.9) 



Comparing (3.9) with the corresponding high-temperature expansion [14] [27] shows that 
up to the order calculated it coincides with AEifi(X) + A£ , 2 ) o(A) at the dual point A = A 0, 
ip with interchanged. In fact, we can also estimate the error we have made in calculating 
(3.9). The error we are making when replacing the true eigenvector by (3.8) is of the 
magnitude iV -1 . This behaviour is preserved when the potential V and the resolvent g 
are applied. Thus, we obtain in all orders starting with the first one a deviation which 
is of magnitude iV -2 . This further supports the identification of (3.9) with the dual of a 
two-particle state. Of course, (3.9) holds only for < ^. 

For > §, the states (3.2) are not the first excited states any more. Now we have to 
consider the following states: 



1 

n 



|| 1 ... 1 2 ... 2 ... ... 0)) P + u Q ||2 ... 2 ... 1 ... 1)) P 
+ . . . + u;^- 1 ) 1 1 . . . . . . (n - 1) . . . (n - 1))) . 
Going through the same steps as before we find: 
lim AEq i =12 sin \^— \ - A4V3cos (^) 

N^oo \ 6 J \3/ 



(3.10) 



2 I cos(^) -2 cos 3) + 1 



7T 



(3.11) coincides with 3AE'i i o(A) at the dual point A = A with interchanged 0, cp in the 
corresponding high-temperature expansion [14] [27] up to the order calculated. 

For = § the states (3.2) and (3.10) are degenerate. However, for large N the dominant 
contribution comes from the states (3.10) such that (3.11) is valid for = ^ as well. At 
= ^ = ^ this is in agreement with the exact result of [10]: 

lim AE Q 1 = 6(1 - A) for0 = <^=^. (3.12) 



4. Duality of spectra 

The results in the previous section remind us of some well known results about duality: 
The spectra in the low-temperature phase at A are dual to those in the high-temperature 
phase at A = A if we interchange and au- This statement for Z n quantum spin chains 
has been known for a long time [28] [29] and was also used in [1]. However, special attention 
has to be paid to the boundary conditions when performing duality transformations. It has 
been observed in [2] that the duality transformation interchanges the role of the charge 
Q and boundary conditions R. In this section we reformulate the precise statement of 
duality for the general Z n -chiral Potts quantum chain (1.1) and discuss its consequences. 
For completeness, a simple non-standard proof is presented in appendix A. 
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We denote the Hamiltonian (1.1) including the parameters by H^(X, R lt , a k f , a^). If 
the Hamiltonoperator is properly normalized for high-temperature expansions (H(X) := 
XH(X~ 1 ) ) it will be called (A, R ht , d 1 ^, «£*), writing again explicitly the corresponding 
parameters. In order to be able to distinguish the parameters we have introduced an upper 
index. Furthermore, abbreviate the space with charge Q ht in the high-temperature phase 
by 1i9 and the eigenspace of f(Q) to eigenvalue by . Now we can formulate the 
statement: H^(X, R ht , restricted to H Qht and (X, R lt , a l £, a^) restricted to 

iiP * have the same spectra if 

Q lt = R ht , R lt = Q ht , 

(4.1) 

rv lt - r» ht J 1 - cv ht \ — \ 

(Xfc — CXfc , (Xfc — CXfc , A — A. 

The momentum decomposition can be applied alike in the high- and low-temperature 
phase. Thus, the statement of duality is also valid if we further restrict to eigenspaces 
with momentum P. 

Note that the duality (4.1) preserves the condition (1.3). Thus, each integrable chiral 
Potts model is dual to exactly one integrable chiral Potts model. The integrable model 
is self-dual iff A = 1 and = ip G ^Z. The point X=l, <fi = (p = exhibits conformal 
invariance [18]. However, this is not true for the other self-dual points. For example the 
point X=l, (f) = (p=^ is not conformally invariant although in its vicinity non-diagonal 
conformal field theories have been used to derive correlation functions [30]. 

From (4.1) we conclude that the quasi-particle interpretation for the high-temperature 
phase of the general Z n -chiral Potts quantum chain [27] can be pulled over to the low- 
temperature phase. The duality transformation interchanges charge sector Q and boundary 
conditions R. Thus, the ground state of the high-temperature phase is mapped to periodic 
boundary conditions R = in the low-temperature phase. However, the fundamental ex- 
citations are mapped to different boundary conditions corresponding to the charge sectors 
Re {1, ... n — 1}. Therefore we observed only composite particle states in section 3. 

It is well known that in the high-temperature phase the limit iV — > oo of the energy 
eigenvalues is independent of the boundary conditions whereas the charge is substantial. 
The duality (4.1) implies that in the low-temperature phase this limit does not depend 
on the charge, but the spectra are clearly different for different boundary conditions. The 
independence from the charge has already been observed in section 3 for some eigenvalues. 
One can also argue directly that this degeneracy holds for the complete spectra, at least 
in the range where the perturbation expansion converges. 

In the perturbation expansion all orders can be organized with respect to the energy 
eigenvalues at zero temperature. All energy eigenspaces - apart from the ground state - 
have a dimension that grows at least with N. We have seen in section 3 for some examples 
that at a fixed order k of the perturbation expansion only a finite number Ne,u of matrix 
elements of the potential V k in the eigenspace to energy E depend on the charge Q. 
This holds generally as we can see from the proof of duality presented in appendix A (in 
particular (A. 2) and (A. 3) ). Thus, the term for energy E at order k in the perturbation 
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expansion has a Q-dependent term that is at most of order (^jf^) 2 as N — > oo. This 
implies that the differences between the charge sectors converge to zero in the large chain 
limit. 

It remains to check that this is also true for the ground state. The contribution of the 
ground state to a perturbation expansion for any energy level other than for the ground 
state energy itself is neglegible. The kth order of the ground state energy per site is 
independent of the charge sector Q if k < N. This is easy to see because in the perturbation 
expansion application of k powers of the potential V to the ground state is projected back 
onto the ground state. Thus, all excitations that are created must be annihilated again, 
or an excited state is proportional to the ground state. The second case is only possible if 
an excitation is carried around the boundary of the chain implying k > N. On the other 
hand, the first possibility yields results that are clearly independent of the charge sector. 
Finally, convergence of the perturbation series for eo implies that the limit iV — > oo is 
independent of the charge sector Q. 

In summary, in the low-temperature regime all charge sectors are degenerate for iV — > oo, 
at least if the perturbation expansion converges. 

5. Correlation functions in the low-temperature regime 

In this section we will apply methods explained in more detail in [27] to the correlation 
functions in the low-temperature phase of the chiral Z3 Potts quantum chain. Note that the 
duality argument of section 4 applies only to the Hamiltonian and not to other operators. 
Thus, quantities like e.g. correlation lengths may be different in these two phases. In fact, 
current knowledge about the correlation functions is restricted to the conjectures presented 
in [22] and the argumentation in favour of a non- vanishing wave vector presented in [16] [17]. 

We study the correlation functions 



(v\ 


r+r 


v) 


(v\ 


v) 



(v\ 


Vx°0 


v) (v 




v) (v 




v) 


(v. 


v) 




(v\ 


v) 2 



C r (x):= ^ * . U| " / , C a (x):= ^ - ^ " ■ (5 _ 1} 



using a perturbative expansion [25] for the ground state | v) from the state | GS; Q) (3.1). 
The expansion of the ground state in powers of A leads to an expansion of the correlation 
functions in powers of A: 



Cr(x) =^A fc C^(x) , C a (x) =^A fc C( fc )(x). (5.2) 



K k 4 k \ 

k=0 k=0 



In the definition (5.1) it is legitimate to omit the one point functions for the operator Y 
because they are zero due to charge conservation. 

Below, we will first give the final results for general angles <p, if. It turns out that the 
result is too complicated to infer the general form of the correlation functions (5.1). Thus, 
we then specialize to the superintegrable case = <p = § and calculate even higher orders. 
By looking for a good fit we try to guess the structure of the correlation functions. With 
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this experience we turn back to the general case and discuss how the correlation functions 
should change for general 0, <p. 

In order to save space we present only the final results for the correlation functions. 
For Cr{x) one obtains, using again the abbreviations C = cos TZ = 1 — AC 2 : 



Cp\x) =1 , 



,(2) 



CT(*) 

1 



6C 2 



{4,o - 1} , 



C$\x) = 0, 

( 3 ) COS^ 



27C 2 



(1 - 4,o) 



2(1 - 16C 2 ) 1 \ 
3ft 2 + 16C 2 / 



+ 5. 



'l + 2C 2 -3i sin(^) 5 



3ft 2 



16C 2 



The first orders of the correlation function C a (x) read as follows: 

Cl°\x)=5 x ,o , Ci 1 \x) = 0, 

' 2 




ci a) (*) = -j{|? 

/ 7 



+ 4,i 



ft + 2C 



W3)/-N_ cosy? [4,0 4,i 
C^J' °- W " 9C \ 2C 2 + ft 



81C 2 



4,0 



16C 2 + K 



;)+«...( 



15 29 
+ 



16C 4 



i6C 2 12ft ft 3 



IE, 2 



1 7 + 44C 2 
"8C 2 + 12ft 2 



(5.4) 

Note that the correlation functions (5.3) and (5.4) do not depend on the charge sector. 

(5.3) suggests that the correlation function Cr(x) tends to a non-zero constant for large 
distances x - in contrast to the correlation functions in the high-temperature phase [27] 
so that that the low-temperature phase is ordered over long ranges. However, beyond this 
general conclusion, it is difficult to guess from (5.3) or (5.4) what might be the behaviour 
even for small x. Thus, we set <fi = ip = calculate four further orders and obtain 



C 



a 



c 



a 



(x) =1 , 
2 

(*0 =g (4,0 - 1) , 



C£\x) = C^\x) = C^\x) = C^\x) = , 
, 7 5 — iV3 

W = £j" (4,0 - 1) + 4,1 , 

0*0 =^T7 i 336 (4,o - 1) + 1604,1 + 605 x , 2 } - {265 Xjl + 205 x , 2 } 



(5.5) 



(x) 



6561 
1 



354294 



354294 



6561 

{12600 (4,o - 1) + 68525 X)1 + 35215 x , 2 + 12254, 3 } 
{960<5 X)1 + 995<5 X , 2 + 5255 X)3 } 
10 



and 



Ci°\x) =5 X , , 



Ci'Hx) = Cf\x) = Cf\x) = Ci 7 Hx) = , 



{-41d x ,o + 144,i + 85 x , 2 } 



(5.6) 



C - 6) ( x ) =77^ {"586^,0 + 147^,! + 1265 x , 2 + 805 x , 3 } , 

^ 8)(X) = 53li41 ^" 9927 ^>° + 2130 ^-i + 1721^,2 + k5 x , 3 + 9105 X)4 } . 

(8) 

Unfortunately, determination of the constant k in exceeded the numerical range of 
our special purpose computer algebra system. 

Up to the order calculated, C a {x) is real for all values of the parameters 0, 9?, A. Cr(x), in 
contrast, has a non- vanishing imaginary part. By analogy to the high-temperature regime 
[27] and from the results in [16] [17] one might expect that Cr(x) is oscillating. Indeed, the 
correlation functions in the superintegrable case <p = V? = § can nicely be fitted by 



C a {x) = ad xfi + be £ , 
C r (x) = m 2 +pe~^ + ^ x . 



(5.7a) 
(5.7b) 



For A G {j, 5, 1 } good fits to (5.5) and (5.6) using (5.7) are given by the values in the 
following table: 



A 




a 


b 




m 2 


L 




P 


0.25 


0.26(1) 


0.89(1) 


0.10(1) 


0.24(3) 


0.9857605 


30 ± 18 


0.26(3) 


0.011(3) 


0.50 


0.41(3) 


0.85(1) 


0.11(1) 


0.37(5) 


0.9381 


30 ± 15 


0.27(1) 


0.05(1) 


0.75 


0.59(9) 


0.75(2) 


0.14(2) 


0.50(7) 


0.832 


30 ± 13 


0.27(1) 


0.15(2) 



Table 1: Parameters for the correlation functions (5.7) at = ip = ^ 

First, we remark that the correlation lengths satisfy £ CT = £r = : £ for all values of A within 
the numerical accuracy. In fact, one expects this equality because the correlation lengths 
should be the inverses of some mass scale, and there is only one mass scale in our problem 
because all three charge sectors are degenerate. Furthermore, we observe that our data 
is compatible with an oscillating correlation function for the operator T. The oscillation 
length L (or wave vector) is around 30 sites in a major part of the low-temperature phase. 
In [16] [17] it has been predicted that L should diverge as A crosses the phase boundary 
and approaches A = 1 where the critical exponent is expected to equal |. Our results are 
compatible with a divergent oscillation length at A = 1 although due to the large errors 
we do not even see that L increases with A. 

We should mention that the linear approximation e~^ ±x rs 1 — ixsin (^) is as good as 
(5.7b). However, (5.7b) seems to be a more natural form. The fact that the relative error 
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of the estimate for q is much smaller than that of L just comes from the fact that L and 
q are related by exponentiation. 

Note that the short correlation lengths strongly damp the correlation functions: At x = 7 
where we expect the first zero of Cr(x) it has already decreased by at least six orders of 
magnitude for A < |. Closer to A = 1 the correlation length should increase but so should 
the oscillation length as well. Thus, it will be very difficult to obtain more precise results 
from approximative arguments and an exact expression for Cr(x) is probably needed in 
order to decide whether (5.7b) really is the correct form and to determine the wave vector 
L accurately. 

Before we conclude the discussion of the correlation functions for the superintegrable chiral 
Potts model, we mention that a conjecture for the form of Cr(x) has been formulated in 

n — x 

[22]: Cr(x) = m +0(e ? r ) where m is the order parameter. Our result (5.7) is compatible 
with this from. In [1] the conjecture for the order parameter 

has been formulated, but (5.8) has not been proven yet. The constant term in (5.5) is in 
exact agreement with (5.7) and (5.8) up to the order calculated, such that we may assume 
that at least the constant term of Cr(x) is now known exactly. 

A few remarks on the choice of ground state (3.1) are in place because in [1][9] (5.8) 
has actually been derived considering an expectation value of the operator r x . We have 
already pointed out that the one point functions of T x vanish identically due to charge 
conservation if one uses the charge eigenstates (3.1). However, if one uses instead non- 
charge eigenstates like | ... 0) for a perturbative expansion of | v) they do not vanish. 
Indeed, using an expansion for | v) from | ... 0) we once again verified equality of this 
one point function with the order parameter m. If we redefine Cr{%) by replacing | v) by 
| v) and subtracting the contribution from the one point functions, this is in fact the only 
change, i.e. Cr(x) = C-p(x) — m 2 . C a {x) remains unchanged under this redefinition. 

For more general values of the angles 0, p one expects the correlation functions to 
be also of the form (5.7) - of course with different values of the parameters. We can see 
from the constant term m 2 of the correlation function Cr(x) (5.5) that it will not be of 
the form (5.8) for general 7^ ^ 7^ cp. In general, the coefficient of A 3 for the constant 
term does not vanish and m 2 does not even have an expansion in powers of A 2 . Among 
the powers that we have calculated for the general case only the fourth order in (5.5) has 
a no n- vanishing imaginary part at x = 1. Under the assumption that (5.7b) is the general 
form we would expect the imaginary part at x = 1 to be proportional to sin (^) for very 
small temperatures A. Thus, we expect for very low temperatures A the relation L _1 ~ 0. 
On the one hand, this explains the conjectured presence of a second length scale L in 
addition to the correlation length £. The oscillation length L just comes from the chiral 
angle and thus these two scales must be related to each other. On the other hand, the 
oscillation (should it really be present) will vanish smoothly as the parity conserving Potts 
case = (p = is approached. 
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6. Conclusion 



In this paper we discussed the low-temperature phase of the ^3-chiral Potts quantum chain. 
We have perturbatively calculated the ground state energy and the first gaps for P = at 
general chiral angles 0, 9?. We explicitly observed duality to the high-temperature phase 
and independence of the charge sector Q. This demonstrates a general duality property 
stating equality of spectra in the low- and high-temperature phase. Thus, a quasi-particle 
interpretation for the high-temperature phase of the general 2 n -chiral Potts quantum chain 
[27] can be pulled over to the low-temperature phase. However, charge Q and boundary 
conditions R are interchanged by the duality transformation. In particular, for periodic 
boundary conditions one sees only energy levels above the ground state that correspond 
to composite particle states. We also gave a general argument that in the infinite chain 
length limit all charge sectors are degenerate. 

We have further studied correlation functions for the operators a and V in the low- 
temperature phase of the 23-chiral Potts quantum chain. The correlation function C-p(x) 
has a constant term m 2 indicating long range order. Fitting (complex) exponential func- 
tions to the perturbation expansions of the correlation functions we estimated the corre- 
lation length £ for = ip = | and found agreement with the prediction [16] [17] that the 
correlation function C-p oscillates. A rough estimate for the oscillation length L has also 
been obtained. We argued that the oscillation length should satisfy L ~ <p~ l for small 
temperatures. In particular, the oscillation vanishes for vanishing chiral angle (f). These 
first results make an exact determination of the correlation functions desirable. 
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Appendix A: A Proof of duality 



Duality has been proved in e.g. [28] [29] . Still, we would like to present a slightly different 
approach in this appendix. We derive duality by comparing the representation f (2.3) to 
the following representation r that is usually considered: 

r(aj) \ h...ij...i N ) = u %i | h . . . % i . . . i N ) 

r(Tj) | %\ . . . ij . . . in) =| i\ . . . (ij + 1 mod n) . . . zjv) • 

Note that r(Tj) = r(aj) and r(aj) = r(T~j~) and that the representations f and r are 
unitarily equivalent. 

We recall the statement of section 4 before proving it. Let Hffi (A, R ht , a 1 ^, a 1 ^) be 
the Hamiltonian with suitable normalization for high-temperature expansions (H(X) '■= 
Ai7(A _1 ) ) and (A, R lt , a%, a l £) be the Hamiltonian (1.1) with corresponding parame- 
ters. Furthermore, abbreviate the eigenspace of r(Q) to eigenvalue by H9 and that 
of f(Q) to eigenvalue u qU by H qU . Then (A, R ht , af , af ) restricted to H qM and 

''(A, R lt , a^, a 1 ^) restricted to H® 1 * have the same spectra if 



Q lt = R ht , R lt = Q ht , 

(4.1) 

„lt _ n ht At _ -ht \ - \ V ' 



For the proof we fix the state | GS; Q lt ) to be the ground state (3.1) in . Then the 
following states are a basis for VP * 1 ): 

\Q;i 2 ...i N ) :=f(ap 2 )...f(a N lN )\GS;Q lt }. (A.2) 

Note that this implies: 

f(ai) \ Q;i 2 ...i N ) = u QU | Q; (i 2 + 1) . . . (i N + 1)) . (A.3) 

Now consider the following intertwining isomorphism /: 

I \Q;i 2 ...i N ) :=\(-i 2 ){i2-i 3 )---(iN-i-iN){iN + R lt ))- (A.4) 

Note that the map (A.4) maps the ground states in both phases onto each other. It is now 
straightforward to check using the basis (A.2) that 

Ir(°(j+i mod n) ) = r(r,T+ +1 )J , 7f(r,T+ +1 ) = r{a 3 )I . (A.5) 

The observation that / is a unitary map and r and f are unitarily equivalent in conjunction 
with (A.5) proves duality. 



l ) these may be different from e.g. (3.2) and (3.10) 
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